The screening of 4/ moments and delocalization in the compressed hght rare earths 
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Spin and charge susceptibilities and the 4/", 4/"~^, and if"'^^ configuration weights are calcu- 
lated for compressed Ce {n— 1), Pr {n — 2), and Nd (n = 3) metals using dynamical mean field theory 
combined with the local-density approximation. At ambient and larger volumes these trivalent rare 
earths are pinned at sharp 4/" configurations, their 4/ moments assume atomic-limiting values, are 
unscreened, and the 4/ charge fiuctuations are small indicating little / state density near the Fermi 
level. Under compresssion there is dramatic screening of the moments and an associated increase in 
both the 4/ charge fluctuations and static charge susceptibility. These changes are coincident with 
growing weights of the 4/"~^ configurations, which it is argued are better measures of delocaliza- 
tion than the 4/"*^ weights which are compromised by an increase in the number of 4/ electrons 
caused by rising 6s, 6p bands. This process is continuous and prolonged as a function of volume, 
with strikingly similarity among the three rare earths, aside from the effects moderating and shift- 
ing to smaller volumes for the heavier members. The observed a-7 collapse in Ce occurs over the 
large-volume half of this evolution, the Pr analog at smaller volumes, and Nd has no collapse. 
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I. INTRODUCTION 

The trivalent rare earth series is an important but 
not well understood test ground for the study of strong 
electron correlation, and the manner in which its ef- 
fects diminish as the 4/ electrons dclocalize under pres- 
sure. These metals appear initially to remain local- 
ized as they pass through a sequence of high-symmetry 
close packed phases keyed to 3d-band occupancy^* and 
then on further compression, eventually reach low- 
symmetry early-actinide-like structures suggestive of /- 
electron bondingi^iii^ Transitions in the region between 
these two limits may exhibit unusually large volume 
changes. While not the case for Nd^^i^ the 7-a "volume 
collapse" of 15% in Ce is well known ^i^iSiiiS and sim- 
ilar collapses occur in Pr (9%)^ii^ii^aiii^ Gd (5%)^ 
and Dy (6%) Magnetic properties at atmospheric pres- 
sure are generally consistent with atomic 4/ Hund's 
rules moments while the susceptibility for the col- 
lapsed a-Ce phase- and for early actinide analogs^^^ is 
temperature- independent, enhanced Pauli paramagnetic, 
indicating absent or screened moments. On the other 
hand, high energy neutron scattering measurements for 
Cer^ and x-ray emission spectroscopy for Gd^^ continue 
to detect 4/ moments in the collapsed phases, possibly 
sensing underlying "bare" moments in spite of screening 
effects. The 4/ electron delocalization itself may be ex- 
amined using resonant inelastic x-ray scattering determi- 
nation of the probabilities of finding /"^^ configurations 
in a compressed rare earth of nominal /" character 

There are at present two viable explanations for the Ce 
collapse, with possible implications for the other triva- 
lent rare earths. One is that it is driven by a Mott 
Transition (MT) in the 4/ electrons, ^^'^^ while the other 
Kondo Volume Collapse (KVC) model points to rapid 
volume-dependent changes in screening of the 4/ mo- 
ments by the valence electrons j^^'^^i^^'^^ The conflict be- 



tween these scenarios is exaggerated by the use of in- 
compatible approximations. Polarized local-density ap- 
proximation (LDA), LDA-I-U, and self-interaction cor- 
rected LDA calculations have been used to support the 
MT picture,^i^i^i^i23i3ii35 While valuable, these are 
still static mean-field treatments which yield either com- 
pletely itinerant (no 4/ Hubbard structure) or completely 
localized (Hubbard splitting but no Fermi-level 4/ struc- 
ture) solutions, thus indicating a too abrupt picture of 
the collapse transitions. On the other hand, the Ander- 
son impurity model treatments^'''' —^'^ used to elucidate 
the KVC scenario can be faulted for omitting direct /- 
/ hybridization and Kondo lattice effects, and one may 
worry whether 0{1/N) solutions^ might favor the local- 
ized limit. 

The combination (LDA+DMFT)^i^ of LDA input 
with truly correlated Dynamical Mean Field Theory 
(DMFT)2£iiS solutions has offered a new perspective 
which has generally been supportive of the KVC sce- 
nario for Ce4ii^i^i^i^!^^ Such calculations for Ce also 
point to ongoing 4/ delocalization in the relevant volume 
range, a critical driver of Mott transitions42ii^ To fur- 
ther clarify the behavior of the compressed trivalent rare 
earths, the present paper reports LDA-fDMFT calcula- 
tions of the 4/ spin and charge susceptibilities and the 
4/", 4/"^^ configuration weights for the first three mem- 
bers, Ce (n = l), Pr (n = 2), and Nd (n = 3). This work 
follows an earlier effort which examined the equation of 
state and spectra for the same materialsi^ Here we con- 
firm that Ce, Pr, and Nd remain localized at pressures 
up through the face centered cubic (fee, 7 for Ce) phases 
as indicated by sharp 4/" populations, unscreened mo- 
ments with atomic-limiting values, and small charge fluc- 
tuations indicating little 4/ state density overlapping the 
Fermi level. On subsequent compression there is rapid 
and dramatic screening of the moments and concurrent 
increase in charge fluctuations and the static charge sus- 
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ceptibility. These changes are also coincident with rapid 
growth in the 4,f"~^ configuration weights, which we ar- 
gue offer a truer measure of delocalization than do the 
4/"+i weights which are complicated by the overall in- 
crease in the number of 4/ electrons due to rising 6s, 6p 
bands. These trends are continuous and prolonged as a 
function of compression, and strikingly similar among the 
three rare earths, suggesting a robust underlying progres- 
sion which must first be acknowledged before tackling in 
general the location or absence of volume collapse transi- 
tions at various stages along the course of this evolution. 

In the remainder of this paper, the susceptibility and 
configuration weight formalisms are first reviewed in 
Sees. II and III, respectively. Computational details are 
given in Sec, IV, results in Sec. V, and a summary in 
Sec. VI. The Appendix discusses the optimal disposition 
of diagonal, one-body /-/ terms used here in the Quan- 
tum Monte Carlo solution of auxiliary impurity problem. 

II. SUSCEPTIBILITY 

The important local susceptibilities for real, multiband 
systems would appear to be associated with the total spin 
S, orbital angular momentum L, total angular momen- 
tum J = L + S, and charge 

Xsir) = (r.S(r).S(O)) (1) 

XL(r) = (r,L(r).L(0)) (2) 

Xj{r) = (T,J(r).J(0)) (3) 

Xc(t) = {Tr[ni{r)~ni][ni{0)-ni]), (4) 

for the electrons in some I shell (e.g., 4/) on a particular 
site, and with T-r the imaginary time r ordering operator. 
In the last, fii is the total number operator J2m a ^rnr^ 
with m — —I, + ,Z, and ni is its generally non- 
integral expectation. For cubic symmetry, 

Xs{t) = 5{TrS,{T) SM) 

= ^{TAMr)-niir)][n^iO)-nim) (5) 

where fia- = ''^■mo- is the total number operator for a 
given spin a = 2 nis = ±1 summed over orbitals m = mi. 
Aside from the factor of 3/4 this is the spin susceptibility 
of the two-band model of Koga et al.,^'^ or for r = 0, 
the same factor times the bare local moment of one- 
band Hubbard and Anderson models. Similarly xl(t) = 
3{TrLz{T)Lz{0)) for cubic symmetry and thus 

Xl(t)==3 ^ mm' {TrnmcriT)n„r>cr'iO)) (6) 

which is 3/4 times the orbital susceptibility of the two- 
band model of Koga et al.j'^'^- taking / = 1/2 and m — 
-1/2, 1/2. 



Spin susceptibility xji''') 

In the presence of the spin-orbit interaction, it may be 
more useful to work in a relativistic basis j ~ I ±1/2 
(except just j — 1/2 for 1 = 0), with magnetic quantum 
numbers v = — j, — j + 1, • • ■ , J. Again for cubic symme- 
try, 

Xj(t)=3 ^^'(r.n,.(T)ri,v,(0)). (7) 

For DMFT calculations which include spin orbit and 
the Hubbard repulsion U, but omit the Hund's rule in- 
traatomic exchange terms, a reasonable approximation 
to the self energy is Yij^^ji y{iuj) ~ 5jji5u_,jiYij{iLo). A 
consistent approximation to {Trhj^{T)fiji,ji{Q)) is 

{Trhj^{T)nj,^>{Q)) 

[ 7V,(T)/(2j + l) if j=/, v = 

~ <^ D,,(t)/[j{2j + 1)] if j=/, (8) 

[ Z),y(T)/[(2j + l)(2/ + l)] ifjV/, 

where 

N,{t) ^ Y.<Trn,,{T)n,M) (9) 

V 

+ n,v(^)"j.(0)]) (10) 

with v ranging over the 2j-|-l states of j (similarly v' and 
j'). Then 

Xj{r) = 5](j- + 1)[jW,(t) - D,,{t)] . (11) 

3 

Note that at t = 0, Nj{0) = nj and Djj'{0)=djj> where 
the number of electrons in the Z-shell is ni = ni+n2 and 
the associated double occupation is di = dii+di2+d22, 
with subscripts 1 and 2 labelling j = 1—1/2 and j = 1+1/2 
for 1>1, respectively. 

The bare or instantaneous local moments correspond- 
ing to each of Eqs. ((lH3]) are given by their t = values, 
e.g., 

MJb + l) = (P)=xj(r = 0) 

^ J2(^ + mjnj-d,,]. (12) 
j 

For / electrons with integer nj = n shell populations 
< n < 14, one might expect in the strongly localized, 
atomic limit that 

ni = min(n,6) 
dii = ni(ni-l)/2 

712 = max(0, n — 6) 
d22 = 712(^2 -l)/2, 
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FIG. 1: Atomic moments (J^) — Ja{Ja+l) for the rare earths. 



which leads to the filled data points (only spin orbit) 
in Fig. 1. Inclusion of the appropriate intraatomic ex- 
change terms would give the correct Hund's rules val- 
ues designated by the open squares. The spin-orbit-only 
results are seen to give the correct qualitative behavior 
with filling, and the correct values of (J^) for subshells 
with one or no holes or electrons. As will be seen in 
this paper, they also appear to give the qualitatively cor- 
rect evolution from localized to itinerant behavior with 
compression, since much of that originates from volume- 
dependent changes in the double occupation which is cap- 
tured correctly. 

Information about screened moments comes from the 
static susceptibility xj{<^ = 0)- Given Curie- Weiss be- 
havior, an effective moment can be extracted from the 
slope of x,7(w = 0) versus T"^, thus 



JsiJs + l)=Txj{LO-- 



(13) 



At large volume and thus weak hybridization, J approxi- 
mately commutes with the Hamiltonian so that Xj(t) ~ 
constant, and thus there is no screening (Js ~ Jf,). There 
is also no screening in the high temperature limit since 
Xj{t) Xj(0) as < r < l/T ^ 0. At small volume 
(strong hybrization) and low temperature, Xji^) istWs 
away between its maximal values at r = and /?, leading 
to Js < Jb- See Fig.H 




FIG. 2: X •'{''') for Ce at 632 K and various volumes. The 
T = values (circles) give the bare moment squared (J^) = 
JbiJb + 1), while the r average gives the screened quantity 
Js(Js-l-l). There is little screening at large volume, while 
the two moments are quite different at small volume. The 
a and 7 sides of the collapse are at volumes of 27.8 and 
33.1 A'^ /atom, respectively, correspondingly roughly to the 
middle two curves. The equilibrium volume is 34.4 A^/atom. 



Charge susceptibility Xc(t") 

The exact expression for the local charge susceptibility 
Eq. Q may be written using the definitions Eqs. (|9ll0p 

as 

Xc{t) = ^^■(^) + 2 E ir) - ■ (14) 

3 3<3' 

Following the language of Ref. IH the local charge "fiuc- 
tuations" are 

{Sfif) EE {{ni~nif)=Xc{T^Q) 

= 2di - ni{ni~ I) , (15) 

while the local static charge "susceptibility" is 



Note that Tx'P < {5n^) since Txc''' is the r average of 
Xc{t) which drops from its t = 0, /? maxima of {5hf) to 
smaller values in the mid t range, similar to Fig. [2] for 
the spin case. This upper bound for Tx^c^ is of interest 
since it helps to identify a large x^c \ which signals the 



drxdr) . 



(1) 



(16) 
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existence of prominent low-energy charge excitations as 
occurs, e.g., in the Yb valence transitions^ 

Clearly {6nf) and Tx'c'^ are the charge susceptibility 
analogs of Jf,(Jb + l) and Js{Js + l), respectively, for the 

spin case. Similarly, (Snf) and Txi^'' must approach one 
another in the large volume localized limit as fii becomes 
an eigenoperator of the system with vanishing hybridiza- 
tion. However, in contrast to the spin case, {Snj) — > 

(and thus also Xc^' ^ 0) in this limit for the trivalent 
rare earths since Uf ^ n and df n(ri — 1)/2, where n 
is the nominal integer 4/ occupation, e.g., n— \ for Ce. 

III. CONFIGURATION WEIGHTS 

The probabilities or configuration weights W]^ of find- 
ing integer k /-shell electrons on a given site are useful in 
discussing delocalization, and are related to Xc{t = ^) = 
{Sfif) insofar as they may also be expressed as linear com- 
binations of ni and di near the localized limit. The Wk 
are given by 

Wk = Zk/^Zk> 

k' 

Zk = ^(fcsfc|e-'^(^-^^)|fcsfe). (17) 

Sk 

Here H is the Hamiltonian; iV, the total number oper- 
ator for all types of electrons; /i, the chemical poten- 
tial; and {|fcsfc)}, a complete set of eigenstates of fii, 
ni\ksk) — k\ksk), where all other quantum numbers be- 
sides k are lumped into Sfc. Evaluating the thermal ex- 
pectations (• ■ • ) of 1, fii, and n;(fi; — 1)/2 using the same 
complete basis yields 

1 = y^Wfc 

k 

ni = fc Wfc 

k 

di = ^fc(fc-l)u;fc/2, (18) 

k 

which shows the statistical nature of ni and di. 

At sufficiently large volumes and low temperatures 
where only Wk for k — n, nil are non-negligible, these 
three Wk may be expressed via Eq. (|18p in terms of ni 
and di 

Wn-i = di-df''{ni) + {\ni-n\-ni+n)/2 

Wn = l-2[di-df''{ni)]-\ni-n\ 
u;„+i = di-df''{ni) + i\ni-n\+ni-n)/2. (19) 

Here it is convenient to use a function d™™(n;) which is 
the minimum possible double occupation for an ensemble 
of sites whose average /-shell population is nj. This is a 
piecewise linear function which assumes the values fc(fc— 
l)/2 at integer k values of ni, and may be expressed. 

dfim) = (n(n-l) + (2n-l)(n;-n) + \ni~n\)/2 (20) 



for the range n — 1 < ni < n + 1. 

Equation appears intuitively to separate the ef- 
fects of delocalization from those arising more simply out 
of changes in ni due to a possible /-shell electron reser- 
voir. Should ni increase due to such a reservoir while the 
system is still in the strongly localized limit, where pre- 
sumably di — (if""(n;), then Wn-i — 0, Wn — l — ni + n, 
Wn+i = ni — n. This suggests Wn~i is untainted by 
such reservoir effects for ni > n, or more generally from 

Eq. nil) 

d,-df\m) = ( ^"+1 ^ ^ " , (21) 

' ^ ^ \ Wn-i n<ni<n + l ' ^ ^ 

which we argue in Sec. V is a useful diagnostic for delo- 
calization. By contrast the local charge fluctuations 

(Shf) = Wn^i + Wn+1 - {ni-nY (22) 

appear to mix delocalization and reservoir effects and so 
are therefore less useful. 



IV. COMPUTATIONAL DETAILS 

The LDA-t-DMFT calculations reported in this paper 
have generally been carried out as in previous work on 
the compressed rare earths j^-''*-''^- All calculations were 
performed for an assumed fee structure, and at a temper- 
ature of 632 K (4 mRy). The spin-orbit interaction was 
included in addition to the scalar part of the 4/ Coulomb 
interaction, i.e., the screened Slater integral = Uf, 
however, not the higher Slater integrals (F'^', /c = 2,4,6) 
which describe the Hund's rules intraatomic exchange. 
While this gives the wrong values for some 4/ moments 
in the localized limit, the volume-dependence accompa- 
nying delocalization of these moments may still be rea- 
sonably captured as this appears to follow from fairly 
general behavior in the evolution of such quantities as 
the double occupancy. As in the earlier papers, the LDA 
contribution to the present work was provided by linear 
muffin-tin orbital calculations in the atomic-sphere ap- 
proximation as described elsewherci^ The auxiliary An- 
derson impurity problem was solved using the Hirsch-Fye 
quantum Monte Carlo (QMC) alogorithm )^^i^° with re- 
sults obtained for L = 80 and 112 time slices extrapolated 
to L = cxD assuming a 1 /L^ dependence. This was unnec- 
essary for x,/(t) where the two L values gave essential 
agreement. The disposition of Ufhf terms between the 
kinetic and interaction parts of the auxiliary impurity 
Hamiltonian in the Hirsch-Fye QMC, and the impact of 
this choice on Trotter corrections, is discussed in the Ap- 
pendix. 

The susceptibilities reported in this work were cal- 
culated within the QMC using Wick's theorem, e.g., 
Eq. (154) of Ref. H^. They were obtained from runs of 
-350, 000 and -100, 000 sweeps for L = 80 and 112, re- 
spectively, using previously converged self-energies to get 
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FIG. 3: (color online). Bare Ji,(J6 + l) = xj(.'^ ~ 0) and 
screened Js(Js + l) = Txjij^ = Q) moments for Ce, Pr, and 
Nd at 632 K, relative to the atomic limiting values Ja{Ja + ^) 
of Fig. [T] for the "only spin orbit" case. Smoothed curves are 
drawn through the points giving L = 80 (circles) and 112 
(triangles) results. 



the input bath Green functions. Error estimates were ob- 
tained from Eq. (5.3) of Ref. [sT] in conjunction with an 
examination of the bin-dependence of the data stored as 
a function of sweep. 



V. RESULTS 

We now present results of the susceptibihty and config- 
uration weight calculations giving first comparisons be- 
tween Ce, Pr, and Nd as a function of volume, then turn- 
ing to insights provided by these results in regard to the 
experimentally observed transitions. Definitions of the 
quantities calculated have been presented in Sees. II and 
III. The volume range studied is from 10 to 50 A'^/atom 
which may be compared to 300-K experimental volumes 
of 14.8, 14.1, and 14.2A3/atom at a pressure of lOOGPa; 
and 34.4, 34.5, and 34.2AVatom at GPa; for Ce^i^ 
Pj,^13j52 g^j^jj Ndf^i^ respectively. 

Figure [3] shows results for J{J+1) corresponding to 
the bare Jb(Jb+l) = Xj{t = ^) and screened Js{Js + ^) = 
Txj(w = 0) moments in Ce, Pr, and Nd, divided by the 
atomic limiting values Ja(>/a + l) = 8.75, 14, and 15.75, 
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FIG. 4: (color online). Local charge fluctuations {5nj) = 

Xc(t = 0) and the static charge susceptibility x;p — Xc(i^ = 0), 
the latter multiplied by T, for Ce, Pr, and Nd at 632 K. Both 
show rapid increases under compression coincident with the 
changes observed in Fig.O Tx)t^ is everywhere more than an 
order of magnitude smaller than ((5n^) reflecting the absence 

of low-energy charge fluctuations, and the ratio Txc^V(<^^/) 
becomes even smaller at small volumes reflecting screening of 
the charge fluctuations. All points are l/L^ extrapolations to 
L = oo, with smoothed curves drawn through them. 



respectively, from Fig. [TJ These are the "only spin or- 
bit" values of that figure, which give the correct J{J+1) 
for Ce, however, are 30 and 36% smaller than the true 
Hund's rules values for Pr and Nd, respectively. It is par- 
ticularly evident for the screened results that the changes 
are most abrupt and occur at the largest volumes for Ce, 
and then successively moderate and shift to smaller vol- 
ume for Pr and then Nd, a pattern which will be seen 
throughout the present results. The bare moments in- 
crease with compression simply because the Uf values 
increasei^ Since a completely random population of the 
4/ states would also have ( J^) and therefore (J^) increase 
with Uf, this does not imply the kind of coherent phys- 
ical moment at the smallest volumes that one certainly 
has in the large volume localized limit. 

Figure 3] shows the local charge fluctuations {Suj) — 
Xc{t = 0) and T times the local static charge suscepti- 
bility Txc''' = Txc{>jJ = 0), which are the charge analogs 
of the bare and screened J{J+1) of Fig. [3] Both quan- 
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titles show rapid increases with compression coincident 
with the changes seen in Fig. [31 Like the spin case 
these two quantities must also approach one another in 
the large volume limit at low temperatures, although 
these limiting values are not finite but for the charge 
case. This follows from nj and df —^ — 1)/2 in 
this limit, where n — 1 (Ce), 2 (Pr), and 3 (Nd), and 
given that (Sfi'j) = 2df — n/(n/ — 1) from Eq. and 

(Snj) > Txc^^ as discussed in Sec. II. The significance of 
the vanishing 4/ charge fluctuations in the large-volume, 
low-temperature limit is disappearance of the Kondo res- 
onance thus leaving a gapped 4/ spectra overlapping the 
Fermi level. 

Also interesting in Fig. 2] is how much smaller the local 
static charge susceptibility Xc"^ is compared to its limit- 
ing maximum {Sfij)/T, e.g., Txi^V(^"/) = 0.017, 0.031, 

and 0.033 at V = 15, 28, and 4lA'^/atom, respectively, 
for Ce, with similarly small values for Pr and Nd. In con- 
trast this ratio is 0.25, 0.58, and 0.09 at the same volumes 
for the second to the last rare earth, Yb^S As pointed 
out in Ref. IH the small ratios for Ce, Pr, and Nd con- 
situte normal behavior where the large onsite Coulomb 
interaction supresses low-energy charge excitations lead- 
ing to a small local static susceptibility Xc^''- The oddball 
is Yb which has a valence transition from divalent (/^^) 
at large volume to trivalent (/^^) at small volume. The 
large Yb ratio 0.58 at 28 A^/atom is in the midst of the 
valence transition where the near degeneracy of /^^ and 
/^^ configurations leads to prominent low-energy charge 
excitations and a consequent large local static charge sus- 
ceptibility. Finally, all four rare earths exhibit decreas- 
ing ratios Txc'^ / (Sn^j:) for compression approaching the 
smallest volumes considered, which reflects screening of 
the charge fluctuations similar to screening of the mo- 
ments as has been noted. 

The probability or configuration weight Wk of finding 
integer k f electrons on a given site is of some interest in 
understanding the manner in which the rare earths evolve 
from localized to itinerant character under compression. 
As all sites are pinned at specific integer occupations n 
in the large-volume, localized limit at low temperatures, 
then for some range of smaller volumes away from this 
limit only Wn and Wn±i are nonneglibile, and may be 
determined from the average number of 4/ electrons nf 
and their associated double occupation df according to 
Eq. p^ . Figure [5] shows Wn-i and Wn+i (the latter +0.1 
for visual clarity) calculated in this manner for Ce, Pr, 
and Nd {n = 1, 2, and 3, respectively). Equation pQ]) 
assumes Wn — 1 — so that the limit w„ = 1, 

Wk^n = is evidently being approached at large volumes. 

As volume is reduced away from the localized limit, two 
things happen. First, for purely one-electron reasons, the 
6s and 6p bands rise relative to the 4/ levels, causing an 
increase in Uf,— and thus a shift in weight from Wn to 
Wriri-i- Second, as hybridization grows, the 4/ electrons 
begin to hop to / states on neighboring sites or into and 
out of valence levels, causing both Wn±i to grow at the 
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FIG. 5: (color online). Configuration weights Wn-i and Wn+i 
(latter +0.1) for Ce (n = l), Pr (n = 2), and Nd (n = 3). The 
former is a measure of delocalization Wn-i = df— d™'°(n/), 
since n < Uf < n + 1 everywhere, and shows a dramatic rise 
with compression for Ce, structure which successively soft- 
ens and shifts to smaller volumes in moving to Pr and then 
Nd. The weights Wn+i are complicated by the general in- 
crease in n/ with compression due to the rising 6s, 6p bands. 
All points are extrapolations to L — oo, with smoothed 

curves drawn through them. 

expense of w„. Only is a true measure of the second, 
delocalization effect uncomplicated by the consequences 
of increasing Uf^. Or more generally, this diagnostic 
for delocalization is given by Eq. (PT|) which is w„_i for 
Uf >n and Wn+i for Uf < n. 

Figure [5] shows clear evidence of both the the ris- 
ing 6s, 6p bands and delocalization. Note first the on- 
set of delocalization as Wn-i increases with compression, 
and how much more abrupt this behavior is compared 
to the smoother Wn^i curves which reflect also the ef- 
fects of increasing n/. These onsets occur near the 
equilibrium volume of Ce (Vb = 34.37 A'^/atom), but at 
somewhat smaller volumes than the corresponding Vq 
of Pr (34.54) and Nd (34.17). At the largest volumes 
(I/>42A3/atom), is an order of magnitude larger 
than w„_i < 0.003 for all three rare earths, suggesting 
the behavior there is predominantly the increase in 
due to the rising 6s, 6p bands. The Fermi level in this 
regime lies midway between the two Hubbard bands, and 
a DMFT calculation with the Hubbard I self-energy — 
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which is incapable of generating the quasiparticle peak — 
has Uf pinned to 1 in the case of Ce for all volumes 
V > 17.6 A'^/atom. This suggests that the large- volume 
shift in weight w„ — > is not mixed valence in the 
sense of the Fermi level moving into the upper Hub- 
bard band, but rather that the rising 6s and 6p bands 
transfer electrons into the 4/ quasiparticle peak, which 
is nonetheless still quite small in this region. 

There have been experimental determinations of the 
Ce configuration weights wo,i,2 using Anderson impurity 
model analyses of photoemission^'' and resonant inelastic 
x-ray scattering2S. data, and a number of LDA+DMFT 
calculationa^'^'^'^i^ of these quantities or related 
rif = 1—W0+W2- All are consistent with wi^l, Wk^i — >0 
in the large volume limit. There is an asymmetry be- 
tween f^-f^ and f^-f^ mixing in solutions of the An- 
derson impurity model which tends to lead to the pre- 
dominant transfer wi — > wo with growing hybridization, 
and thus a decrease in across the 7 to a collapse and 
wq > W2 in the a phasci^ This behavior is seen in the 
analyses of both experimental papers.—iSl The present 
and our earlier work^i^i^ concur with the predominant 
transfer wi — >z«o across the collapse, however, has every- 
where a larger W2 than these analyses, quite possibly re- 
flecting the impact onuf of the rising 6s, Qp bands which 
may not be well treated in the impurity model simula- 
tions. Thus we see a non-monotonic dip in across the 
collapse, but generally f^/ > 1 and wq < W2-^^ 

We turn now to possible insights provided by the 
present results into the experimentally observed volume 
collapse transitions (shaded regions) and phase structure 
(demarcated by vertical solid lines) as shown in Fig. [5] 
for Ce^iSiia Pj.^12j^j14j^ Ndj^ At large volume only 
the face centered cubic (fee) phase is labelled, which is 
the end of the localized, trivalent rare earth series,^ with 
all three rare earths found in the preceeding, double- 
hexagonal close packed phase at ambient conditions. The 
two most dramatic theoretical diagnostics are shown for 
comparison, Js{Js+l) and w„_i, and while they were cal- 
culated everywhere assuming an fee structure, it is hoped 
nonetheless that their volume dependence is sufficiently 
insensitive to structure so as to still provide useful in- 
sights. This is not an issue for Ce where both a and 7 
phases bounding the volume collapse are fee. 

LDA+DMFT calculations^i^>i^iMd5d6 have consis- 
tently supported the KVC scenario ^^i^^ '^'^ for Ce. They 
show, e.g., a rapid build up of the Kondo resonance for 
volume reduced across the 7-a two-phase region, which 
is taken as a signature of the onset of screening. The 
present results now directly report the screened mo- 
ment, and indeed one sees a 60% reduction in Js{Js + l) 
from the 7 to the a side of the collapse. Moreover, if 
the Kondo temperature Tk is defined by the condition 
Tx.iiuJ = 0)/xj(t = 0) = 0.5, then Tk = 632 K at the 
volume F = 29.0 A^/atom (vertical dashed line) for the 
isothermal results in Fig. [S^a). This value is consistent 
with the combined experimental results of Refs. [l^ and 
[27l . except for the middle 10-kbar point of the former 
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FIG. 6: (color online). Screened Js{Js + l) — Txj{i^~0) mo- 
ments (solid curves) and configuration weights Wn-i (dashed 
curves) for (a) Ce, (b) Pr, and (c) Nd at 632 K. The for- 
mer are divided by the atomic limiting values Ja{Ja + l) of 
Fig.[T]for the "only spin orbit" case. The smoothed curves are 
drawn through the combined L = 80 (circles) and 112 (trian- 
gles) points for the moments, and extrapolated to L — > cx3 for 
the Wn-i- All calculations were for an assumed fee structure, 
however, the results are compared to the observed phases de- 
marcated by the vertical lines, with shading for Ce and Pr 
indicating significant two phase regions. The vertical dashed 
lines show the volumes where the Kondo temperature is 632 
K. 



which appears out of place given the roughly exponential 
behavior expected for Tk{V)^ 

The competing MT scenario for the Ce volume col- 
lapse relies on 4/ electron dclocalization as the underly- 
ing driving mechanism. ^■^'^'^ If Eq. pT|) is accepted as one 
suitable diagnostic, then wq in Fig.[SKa) suggests that the 
collapse also coincides with dclocalization, a point that 
has been made previously.^^'"'^ To be careful, note that 
Eq. (|19p for the Wk does break down when configurations 
Z*^ for |fc— rt| > 2 become important as the itinerant limit 
is approached, ft would seem from the present and ear- 
lier LDA-I-DMFT calculations that the following are all 
different facets of the same continuous and extended evo- 
lution with compression of a nominally /" trivalent rare 
earth: (1) transfer of 4/ spectral weight from Hubbard 
side bands to a Fermi-level structure, (2) screening of the 
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4/ moments, (3) growth of 4/ charge fluctuations, and 
(4) dispersal of configuration weight away from Wn — 1 to 
neighboring and then more distant configurations. Some 
of these facets look more Kondo like, and some seem more 
consistent with intuitive ideas of delocalization. 

There is striking similarity between the theoretical di- 
agnostics for Ce and those for Pr and Nd in Fig. [51 al- 
though the screening and delocalization effects shift to 
smaller volume and become more gradual for Pr, and 
then more so for Nd. While the Ce evolution seems to 
have saturated under compression on reaching the body 
centered tegragonal (bet) phase at about F = 22 A'^/atom 
(P= 13 GPa),^ Nd has yet to reach this stage even by the 
smallest volume considered 10 A^/atom (P~365 GPa),™ 
and Pr is intermediate. The stability fields of the dis- 
torted fee (dfcc) structures in both Pr and Nd appear 
analogous to the a-7 two phase region in Ce, according 
to progression of the two theoretical diagnostics. Only 
on further compression is there a volume collapse in Pr 
from dfcc to a-U, while Nd passes through two addi- 
tional phases before reaching the same a-U structure, 
absent any large volume changes. It is conceivable that 
this progression in collapse size and location follows sim- 
ply from the shift to smaller volumes and moderation 
in the correlation contributions which must then com- 
pete with the ever bigger underlying benign part of the 
equation of state. While there are suggestions of a Van 
der Waals loop in the LDA+DMFT free energy corre- 
sponding to the isostructural fee Ce collapse r^''^'^'^ 
one must worry about the need to include the proper 
Hund's rules exchange for multi-/ electron Pr and Nd, 
as well as performing the calculations for all of the ob- 
served structures, e.g. a-U. Moreover, the Ce collapse 
has a critical temperature of only 480 while for Pr, 
the dfcc phase is absent above about TOOK with yet a 
new Pr-VI phase intermediate between the a-U and fee 
phasesj^ These temperature senstivities are also a re- 
minder of the need to include lattice vibrational contri- 
butions, which may themselves further modify the nature 
of the collapse transitionsj ^^'^^ 



VI. SUMMARY 



incident with rapid growth in the 4/"'"-'^ configuration 
weights, which we argue offer a truer measure of delo- 
calization than do the 4/"+^ weights which are compli- 
cated by the overall increase in the number of 4/ elec- 
trons due to the rising 6s, 6p bands. Combined with 
earlier LDA+DMFT results, this work suggests a con- 
tinuous and extended evolution with compression of a 
nominally 4/" trivalent rare earth in which there is (1) 
transfer of 4/ spectral weight from Hubbard side bands 
to the vicinity of the Fermi level, (2) screening of the 
4/ moments, (3) growth of 4/ charge fluctuations, and 
(4) dispersal of configuration weight away from 4/" to 
adjacent and then more distant configurations 4/*^. The 
static charge susceptibility Xc"^ mirrors the volume de- 
pendence of the charge fluctuations {Sn'j), and further 
indicates screening of these fluctuations at small volume, 
and supression of low-energy charge fluctuations by the 
Coulomb interaction at all volumes4^ 

The composite evolution in Ce persists until the bet 
phase (y < 12 A^/atom, F > 13 GPa),* whereas the a-7 
two-phase region coincides with only the large-volume 
half or so of this progression. These effects shift to 
smaller volume, become more gradual, and extend over 
a broader volume range for Pr, and even more so for 
Nd, which has yet to saturate by lOA'^/atom (P ~ 
365 GPa).^ The observed stability fields of the dfcc 
phases in Pr and Nd also coincide with the large- volume 
half of the evolution in these materials, and thus corre- 
spond to the a-7 two-phase region in Ce. The collapse is 
observed in Pr only on further compression, and absent 
in Nd. Since LDA-I-DMFT calculations do hint at a Van 
der Waals loop in the Ce free energy at about the right 
place ^i^i^i^ the full story for Pr and Nd may await 
inclusion of the Hund's rules exchange for these multi-/ 
electron cases, assumption of the correct structures for 
all phases (e.g., a-U), and possibly also adding lattice vi- 
brational contributions!^'^ Even so, this work suggests 
that the collapse transitions may not be the predominant 
story in the compressed trivalent rare earths, but rather 
consequences of an underlying and more robust evolution 
associated with 4/-electron delocalization. 



We have reported LDA+DMFT calculations as a func- 
tion of volume at 632 K, for the 4/ spin and charge sus- 
ceptibilities, and the probabilities of finding 4/"^^ con- 
figurations in the nominally 4/" trivalent rare earths Ce 
(n = 1), Pr (n = 2), and Nd {n = 3). We find these 
metals to remain localized at pressures up through the 
fee (7-Ce) phases, the last structure of the initial close- 
packed series,^ as indicated by sharp 4/" populations, 
unscreened moments with atomic-limiting values, and 
small charge fluctuations indicating little 4/ state density 
overlapping the Fermi level. On subsequent compression 
there is rapid and dramatic screening of the moments 
and concurrent increase in charge fluctuations and the 
static charge susceptibility. These changes are also eo- 



ACKNOWLEDGEMENTS 

Work by AKM was performed under the auspices of the 
U.S. Department of Energy by Lawrence Livermore Na- 
tional Laboratory under Contract DE-AC52-07NA27344. 
Work by MJ and RTS was supported by the SciDAG 
program, grant DOE-DE-FC0206ER25793. This work 
has also benefited from an alliance with members of the 
DOE/BES funded Computational Materials Science Net- 
work Cooperative Research Team on "Predictive Capa- 
bility for Strongly Correlated Systems." AKM is grateful 
for conversations with K. Held, J. Kunes, W. E. Pickett, 
and E. R. Ylvisaker. 



9 



APPENDIX: TROTTER CORRECTIONS AND 
THE BATH GREEN FUNCTION 

One issue related to application of Hirsch-Fye^ QMC 
to the rare earth series deserves attention. To prepare for 
the Hubbard-Stratonovich transformation, one rewrites 
the interaction part / of the auxiliary impurity Hamilto- 



Uf 



a<a' 



Q<a' 



[nafla'-kina+na')] 



(23) 



where a labels the 14 / states. A strict generalization of 
the Hirsch-Fye treatment would be to remove the Uffif 
term from the right side of in Eq. ((23|) and add it to 
the one-body or kinetic energy part K of the Hamil- 
tonian, with the consequent changes in unperturbed or 
bath Green function G(t) and also the Trotter breakup. 
However, for Uf — 6eV and T = 600 K, this would lead 
to an / bath Green function, absent hybridization and 
spin orbit for simplicity, of 



^-Y.Tr[pUr)fi{Q)]/Tr[p] 

-10-^""^/^ (24) 



14 

a 

_g-r(£^6.5C//-p) 



g--/3(e/+6.5(7/-M) + 1 



Here p = exp[—f3{K ~ pN)], p is the chemical potential, 
N the total electron number operator, < t < /3, and we 
take a Ce like site energy e f—p ~ —Uf/2. Such a function 
may get into a region of underflow errors at large r in 
numerical computation. Note this problem arises here 
in the seven-band case only because the Ufhf term in 
question is 13 times larger than in the familiar one-band 
case. For this reason we left the Uffif term alone in 
previous work on the early rare earths r^ii^di and took 
the kinetic and interaction parts of the auxiliary impurity 
Hamiltonian to be 



A, A' 



(25) 



fUfhf , (26) 



Q<a' 



where A ranges over both / and the additional bath de- 
grees of freedom, and Ai covers the rest of the one-body 
terms. Keeping the Uffif term in / can be handled with 
a minor modification of the Hirsch-Fye technique. 

Since Sf—p itself varies from about — 0.5J7/ for Ce to 
about —6.5Uf for Lu, the choice Eq. (|25p then leads to 
similar problems for the bath Green function of the late 
rare earths. A more consistent treatment for the series as 
a whole would be to use one of the Hubbard bands as the 
effective / site energy in the bath Green function, and we 
use the approximate position of the lower Hubbard band, 
Elh = ^/ + (^^l)t^ / j with the integer n being the nominal 
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FIG. 7: (color online). Trotter corrections in (a) xj(''' = 0) = 
(J^) and (b) Xc(''" = 0) — {Snj) for Nd as a function of volume. 
Results for L = 80 and 112 time slices are shown, using the e/ 
(circles, up triangles) and elh (squares, down triangles) effec- 
tive / site energies in the QMC. The lines show extrapo- 
lation to L = oo. It is evident that the Trotter corrections are 
considerably smaller using the elh choice Eqs. (|27I28|) . The 
horizontal dashed line in (a) is the atomic limiting value. 



/ occupation. Thus we take 



K 

i 



[sf + in-l)Uf] hf + J2clAt 



XX' Cy 



A. A' 



a<a' 

+ if-n)Ufnf, 



(27) 



(28) 



which was successfully used in Hirsch-Fye calclations for 
the second to last rare earth, Yb (n = 13),^^ while for Ce 
(n = l), Eqs. (|27I28P are identical to Eqs. (|25I26|) . 

We note here that even for the light rare earths such 
as Pr {n = 2) and Nd {n = 3) that Eqs. (I27I28P serve to 
diminish the Trotter corrections relative to Eqs. (|25l26p . 
This is seen for Nd in Fig. [71 where the instantaneous spin 
Xj(t = 0) = J&(Jfc-fl) and charge Xc{t = 0) = {Sh'j) suscep- 
tibilities are shown as functions of volume. Results using 
Ef [Eqs. (|25l26p . open symbols] and elh [Eqs. (|27l28p . 
closed symbols] as effective / site energies are given for 
L=:80 and 112 time slices. The former points are quite 
different indicating significant Trotter corrections, while 
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the latter are essentially on top of one another except 
for the charge case at the smallest volumes. It is re- 
assuring that extrapolations to L = oo (lines) for 
the Ef and elh cases are fairly close at smaller volumes. 
On the other hand, extrapolation of the e/ results yields 
an unphysical negative Xc-{t — 0) large volumes (not 



shown) indicating non behavior and the need for 
larger L. By contrast the ^lh results are already con- 
verged by L = 80 for both spin and charge cases at large 
volume. It is clear that Eqs. (|27l28p represent the supe- 
rior approach, and all results in the present paper have 
been obtained in this manner. 
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